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We investigate the dc current-voltage characteristics of d-wave Josephson 
junctions, where the barrier at the interface may have arbitrary strength. 
Dividing the current into n-particle currents In (n integer), we can explic- 
itly show which physical processes are responsible for the subharmonic gap 
structure (SGS). For orientations where midgap states (MGS) exist, the res- 
onances in the n-particle processes are drastically changed, giving rise to a 
strongly modified SGS. Introducing broadening in a phenomenological way, 
we show that MGS may produce a current peak near zero bias and we explain 
which physical processes are contributing to this peak. The agreement of our 
theory with recent experiments is discussed. 

PACS numbers: 74.50.+r, 74.25 



1. INTRODUCTION 

The formation of midgap stateJll (MGS) at surfaces and interfaces of 
d-wave superconductors affects the current transport properties of junctions 
involving d-wave superconductors. It has been estabhshedBS that the zero- 
bias conductance peak (ZBCP) seen in normal metal/d-wave superconductor 
(N/d) junctions, are due to the MGS. For the ac Josephson effect it was 
shownBu that subharmonic gap structure (SGS) is in general different in d- 
wave junctions compared to s-wave junctions. In Ref. Q it was shown that 
MGS produce current peaks at voltages of the order of the maximum gap, 
while in Ref. |5| also a singularity at zero bias was found. Later on, in Ref. I, 
it was pointed out that this structure was not seen in Ref. 0, because the 
scattering theory did not include broadening effects. Numerical calculationsQ 
including broadening confirmed the picture outlined in Ref. |^. Despite all the 
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interest in the ac Josephson effect for d-wave junctions, the physics behind 
the current peaks has not been fully understood. 

Recently, a scattering theory based approachB'B has successfully ex- 
plained SGS in conventional (s-wave) junctions for arbitrary transparency of 
the barrier: it agrees for one-channel junctions with experimental^ without 
fitting parameters. In a recent reformulationfiJ of the theory it was shown 
that only currents from multi-particle processes creating real excitations need 
to be summed up, since non-physical currents, present in the original formu- 
lation of the theory, cancels. This proves that the Pauli exclusion principle 
is fulfilled. 

In this paper we extend the reformulated scattering theory described 
above to junctions of d-wave superconductors and provide a deeper analysis 
of the SGS and the effects of MGS. 



2. EXPRESSION FOR THE CURRENT 

Considering transport in the a6-plane, we model the da^Jda^ junction 
[a-L/R is the orientation angle of the left/right superconductort3) as described 
in Ref. 12. This reference also provides a detailed description of the method 
we use to solve the time-dependent Bogoliubov-de Gennes (BdC) equation. 
A quasiparticle incident on the junction at energy E undergoes multiple 
Andreev reflections and builds up a scattering state with amplitudes at the 
sideband energies En = E + neV (V is the voltage; n integer). It can be 
shown that the probability current leaking out at sideband En determines 
the n-particle current /„. At zero temperature we havetJO 

IdciV) = do Vn/„(y), In{V) = decosO f dElP{e,E), 



|aoP)|ao?'o-P 



Q={/,r} 

X [(1 - \an?) + (1 - |a„p)|a„r„+|2] , (1) 

Q ^ ^nO 

(1 - aoaoro-r„o)(l - a„a„r„+f„o) - aoaoa„a„ro_r„+t„oinO ' 
where = ekpLy/lTrh. Above, oq and Un are the Andreev reflection am- 
plitudes for the angle 9 at energies £^0 and £'„ respectively (the barred am- 
plitudes are calculated at 9 = n — 9); ro_ and r„+ describe reflections from 
minus and plus infinity in energy space; tno, tno, rno, and f„o are the ele- 
ments of the scattering matrix describing the region between the injection 
point and the exit point. The weight n of /„ appears because the n-particle 
process involves an effective transfer of the charge ne over the junction. 
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Fig. 1. Contributions from the first four n-particle currents to SGS for the do /do 
junction in (a), and the do/d^^ junction in (b). The dashed lines are the total 
currents. In (b) broadening 77 = 0.01 has been introduced reveahng a current peak 
near zero bias. The angle averaged junction transparency is I? = 0.026 and we 
assume zero temperature. Note the inverse voltage scale in (a). 



3. RESULTS AND DISCUSSION 

Analyzing the current in Eq. ([l|), we see that resonances may appear in 
the propagator GnO) which describes the transmission process from energy 
E to En- There are two types of resonances: bare transmission resonances 
(appearing in t^o) ^-nd boundary resonances (due to the denominator). A 
bare transmission resonance may appear when the trajectory hits a bound 
state at the IS interface, so called de Gennes state. In a short do/do junction 
(as in the short s-wave junction) the de Gennes states are located at the gap 
edges^iving rise to the usual SGS at the voltages eV = 2A/n. In the s-wave 
caseBiy the SGS is due to the combination of onset of the n-particle current 
and resonances in the n + 1-particle current (due to an overlap of a bare 
transmission resonance and a boundary resonance) and the n + 2-particle 
current (two overlapping bare transmission resonances). For the do/ do case 
[see Fig. |l](a)] the physics behind the SGS is the same, but for two reasons 
the structure is smeared and suppressed. First, there are no real onsets 
(the d-wave gap has nodes), meaning that the 1-particle current background 
dominates at all voltages. In addition, the resonances in the higher order 
currents are not sharp because of angular averaging. 

Rotating the right d-wave gap away from the a/j = orientation, the 
de Gennes states are moved from the gap edges to zero energy (MGS) for 
those angles where the gap changes sign after normal reflection at the junc- 
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Fig. 2. Introducing broadening reveals a current peak near zero bias in the ^45/^45 
junction. Decreasing 77, sharpens the peak and moves it to lower voltage. In the 
limit ?7 ^ (dashed line) the peak becomes a delta spike at = 0. In (b) wc show 
that particle currents of high order contribute to the peak. The angle averaged 
junction transparency is D = 0.026 and we assume zero temperature. 



tion. The bare transmission resonance is then moved and, in addition, the 
boundary resonances are lost for quasiparticles injected from the right super- 
conductor. This result in a drastically changed SGS: for the d^jdi^^ junction 
[see Fig. |l|(b)] the only surviving structure is at eV = Aq and it is mainly 
due to the 2-particle current: an overlap between the bare MGS resonance 
and boundary resonances (near the left gap edges) produces the peak. This 
happens also in the 5/^45 junctionll^. We introduce, on a phenomenological 
level, inelastic scattering into the problem by adding a small imaginary part 
ir] to the quasiparticle energy, which results in broadening of all resonances. 
For the do/d4r, junction, a small current peak is then revealed near zero bias, 
as seen in Fig. |l|(b). The peak is due to resonances in particle currents of 
order n = 2 and higher and was therefore not discussed in connection to the 
tunnel limit calculations in Ref. H and |^. 

When MGS are present on both sides of the junction (the ^45/^45 junc- 
tion) boundary resonances can never appear. Consequently current peaks 
are not seen in the IV-characteristics§, leaving only an onset of the current 
at eV = Aq (due to the bare MGS resonance). Again, introducing broaden- 
ing reveals a peak near zero bias, see Fig. |2|(a). In Fig. §(b) we show that 
processes of many orders are, i n this case contributing to the peak. 

In recent experimentsl^llB on bicrystal grain boundary junctions of hole- 
doped cuprates a ZBCP was seen for all orientations of the superconductors. 
No real SGS was seen, only a gap- like structure at eV = Aq. Our predicitions 
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agree with the experimental results, apart from that negative differential 
conductance was not found in the experiments. In another experiment^ 
weak SGS (altough not perfectly at eV = 2Ao/n) was observed for edge 
junctions with the do /do orientation, as we also report here. 

4. SUMMARY 

Dividing the current into n-particle currents, we have shown which 
physical processes are giving rise to SGS and current peaks near zero bias. 
For orientations where no MGS are present in the junction the SGS is at 
eV = 2Ao/n as in the s-wave case. When MGS are present on at least one 
side of the junction SGS is lost and a current peak appears near zero bias 
(if we include broadening into the formalism). Currents of high orders are 
contributing to this peak. 
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